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' ' Abstract 

' By applying analyticity and single channel unitarity we derive a new for- 

■ mula which is useful to analyze the role of the left-hand singularities in hadron 

^ , form factors and in the determination of the resonance parameters. Chiral 

perturbation theory is used to estimate the left-hand cut effects in vrvr scat- 
tering processes. We find that in the IJ = 11 channel the left-hand cut effect 
'nI" , is negligible and in the IJ = 20 channel the phase shift is dominated by the 

\ left-hand cut effect. In the IJ = 00 channel the left-hand cut contribution 

l/^ ' to the phase shift has the wrong sign comparing with the experimental data 

. and therefore it necessitates the a resonance. The new experimental results 

\^ \ from the E865 collaboration is crucial in reducing the uncertainty in the de- 

' termination of the mass and width of the a resonance within our scheme. 

S : PACS numbers: 11.55.Bq, 14.40.Cs, 12.39.Fe 

O \ Key words: dispersion relation, left-hand cut, a meson, chiral perturbation theory 
^ ■ 

p^' 1 Introduction 

(D : 

The dynamical origin of the lightest scalar meson is a subject of long lasting debates 
\ and controversies that were reflected in the name /o(400 — 1200) given to it in latest 

X ; PDG pubhcation |T], §. There is a growing number of studies claiming that the vrvr 

c3 ' scattering phase in the scalar-isoscalar channel J^^ = 0^^ / = which rises steadily 

between the vrvr threshold and the /o(980) region is supported by a broad resonance 
known as the sigma meson. The present situation is demonstrated in Fig.|l] showing 
the position of the a pole in the complex mass plane found in different analysis |^ 
- [|T8|. For most recent studies on related issue, one is referred to Ref. |T9| . 

However, owing to the strong interaction nature, it is difficult to discuss the 
dynamics in the IJ = 00 channel without heavily relying on different models and 
approximations. These models or approximations often violate some fundamental 
principles. For example, in most applications of the K matrix approach the dynami- 
cal effects of the left-hand cut {l.h.c.) are approximated by background polynomials 
which are beyond control theoretically. Strictly speaking, such an approximation vi- 
olates crossing symmetry. It is not always clear how these assumptions and approxi- 
mations influence the conclusion at quantitative or even qualitative level. Therefore 
it is very important to study the problem from different angles and to reduce the 
model dependence as much as possible. 
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Figure 1: The poles of the S"- matrix in the complex mass plane (GeV) corresponding 
to the lightest scalar resonance according to Refs. - [jl8| . 



The goal of this paper is to elucidate the role of the left-hand singularities in 
the determination of the position of the a meson. Our philosophy in this paper is to 
start strictly from first principles and to reduce model dependence as much as pos- 
sible. Throughout the text the main assumptions we use are analyticity and single 
channel unitarity. We will also make use of chiral perturbation theory to estimate 
the left-hand cut contributions to tttt scattering processes in the phenomenological 
application. Chiral perturbation theory encounters problems at high energies, an 
unitarized chiral approach in studying the nn system or even the tttt, KK coupled- 
channel system has been proposed and discussed in recent years with remarkable 
success (see Ref. HOf-l^, and also Ref. 0]). However, our treatment of chiral 
perturbation theory is rather different from the standard applications based on it. 
Essentially we only need its prediction on cuts. We are convinced, by analyzing 
the IJ = 00,20 and 11 channels simultaneously, that chiral perturbation theory 
can be used reliably to elucidate the role of the l.h.c. in the determination of the a 
resonance. 

This paper is organized as follows: In the begining we will establish a new 
representation for a hadron form factor in the single-channel case. We prefer to 
introduce it firstly in the more traditional scattering theory in sec. 0, which may 
benefit some readers. The analytic structure of the vr hadronic form factor and of the 
TTTT scattering S matrix in the context of field theory is analyzed in the new scheme 
in sec. |[ Left-hand cut effects of the tttt scattering amplitudes are estimated in 
sec. ||. The sec. || is devoted to study the pole position of the a resonance. The 
analysis in the IJ = 20 and 11 channel is also made. The sec. 13 is for the conclusion. 



2 



2 The dispersion relation for a scalar form factor 



2.1 Basic formulas 

In this section we consider the relationship between an S'-matrix and a scalar form 
factor in a single-channel scattering problem. For the benefit of the reader, we 
begin with a brief summary of well known properties of the scattering matrix. The 
S'-matrix is related to the partial wave amplitude (we consider the case of the S-wave 
scattering) by 

S{k) = l + 2ikf{k) (1) 

where k is the channel momentum. The one-channel unitarity of the 5* matrix in 
the physical region A; > has the form 

sik)* = siky' , (2) 

Im/(A;) = k\fik)\' . (3) 

The reflection property of the S'-matrix (see e.g. I^Sj) leads to the following relation 
between positive and negative channel momenta k: 

S{-k) = S*{k) (4) 

where the connection between A; < and k > is via an analytical continuation 
in the upper half plane. The Riemann surfaces of the S'-matrix and the scattering 
amplitude as functions of energy E = k"^ have two sheets: 

crM - j Si{E) = l + 2zki{E)fi{E) , lmfc>0 . . 

""^"^ \ SniE) = l + 2tkjj{E)fniE) , Im < 

Here f{E) is the scattering amplitude, the subscript / or // denotes its value on the 
sheet I or II, and ki{E) = —kii{E). In the following, however, we often drop the 
subscript (or superscript) / when it causes no confusion. The reflection property 
d^) together with the unitarity relation (|]) gives the following relations between the 
different branches of S{E) and f{E) 

S„{E) = S',(E) = -gl^^ (6) 

ME) ^ II (7) 

which extend by an analytical continuation from the physical region > to the 
corresponding domains of analyticity on the sheets I and II. 

We define the scalar form factor A{k) as an analytical function of k which satisfies 
the following unitarity relation (compare with Eq.(^) 

Im A{k) = kA{k)f{ky , k>0 (8) 
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and has no singularities, except for possible bound states, in any finite part of the 
upper half plane Im k > (that is the sheet I as the function of energy E). To define 
the scalar form factor uniquely additional constraints are needed. One is a trivial 
normalization condition A{kQ) = Aq with some choice of /cq Aq convenient for 
physical applications. 

Equation (^) has the well known Omnes-Muskhelishvili (OM) solution which 
determines the form factor Aj[E) on the sheet I 



ME) ^ P(EW,[^f;^^iE>^ 



(£>) 



where 5{E) is the scattering phase: 

Si{E) = e'*'^^^) , E>0 (10) 

and P{E) is an arbitrary real polynomial {P{E) is real for real E). Here an once- 
subtraction form of the OM solution is used, and the normalization condition is 
absorbed in the polynomial P{E). 

It is useful to remember that in potential scattering the S'-matrix and the form 
factor can be conveniently defined in terms of the Jost function J-'{k) ||2^ 



Sik) ^ ^ ill) 

m ^ ^ (12) 

where J^{k) is analytical in the upper half plane (Im A; > 0) and has the reflection 
property: 

^(-fc) = J^*{k*) (13) 

In this case the polynomial P{E) is reduced to a constant. Also one gets the relation 
between the values of the scalar form factor on the sheets I and II: 

A.m ^ III ,u) 



2.2 A new representation of the scalar form factor 

In this section we derive another useful representation of the scalar form factor which, 
to the best of our knowledge, was not described in the literature. The starting point 
is the unitarity condition (H) which is written in the form 

ImA,(g)^ -^''^ + "'--^''^-"' = k,(E)ME)ME) (15) 

for E >{). Here the identity fj{E) = fii{E) for real E is used. Since the scalar 
form factor Ai{E) is analytical on the sheet I except for the poles corresponding to 
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Figure 2: The integration contours used in the dispersion relations for the scalar 
form factor. 

bound states at < and the kinematical cut along the real E axis, the following 
once-subtracted dispersion relation can be used: 

ME) -MA) ^ (16) 

^ ' 2m Jo {z-E){z-A) 

(E-A) f kj{z)Aj{z)fn{z)_^^ ^^^^ 



271 Jcr {z-E){z-A) 

where the sum is taken over all bound states with energies are the corre- 

sponding residues of the form factor in the poles on the Sheet I, A is a subtraction 
point, and the contour Cr envelopes the kinematical cut < < oo as shown in 
Fig.0. 

The integration contour can be deformed, assuming that the integrand falls off 
fast enough at the infinity (that is why the subtraction is needed), in such a way 
that it goes around the left hand cut —oo<E<Ec corresponding to the dynamical 
singularities related to the scattering amplitude fii{E) on the sheet II (see Fig.|^). 
In the process of deformation the contour crosses several poles of the integrand. The 
terms {z — E) and {z — A) in the denominator of Eq.(|l7|) correspond to two such 
poles. The bound states, if they exist, contribute via the poles on the sheet I in the 
factor Ai{E') . The virtual and resonant states correspond to the poles on the sheet 
II of scattering amplitude fii{E'). Collecting all these terms we get 

I{E) = Hk,(E)A,{E)f„(E)-k,{A)A,(A)f„{A) + B{E) + R{E)) (18) 
(E-A) f k,(.)A,(z)!,.{A^^ _ j^gj 



2tt Jcl {z-E){z-A) 
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''^''^ = ? {E,-E)iE,-A) ' (20) 

^(^^ = 4^ (i5.-i5)(i5.-A) • (''^ 

In above the function B{E) is given by the sum over all bound states (the poles Eb 
of the scattering amplitude on the sheet I), and the function R{E) is given by the 
sum over all resonances and virtual states (the poles E^ of the scattering amplitude 
on the sheet II). According to Eq.(|^ the poles of Sj{E) correspond to the zeroes of 
Sii{E), therefore for all bound states b 

kj{Eb)fn{Eb) = (22) 
The residues of fn{E) in the poles Ei on the sheet II are calculated with Eqs.(0|l^): 

k,{EMAEr)Res[fn{E^)] = '^^^ = (23) 

Pr = ResAniEr) (24) 
Using Eqs. ([l6| - p4D we get the following representation of the scalar form factor: 

(^-^)^" ^(^) (25) 



TT J-oo {Z - E){Z - A) 

where the function A{E) is defined on the left hand cut — oo < E < El (see Fig.^ 
by 

^ kIiE)AJiE)ifniE + ^e)-fn{E-^e)) ^^6) 

AnjE + ze) - An{E - ze) 

2z ^ 

Here we used Eq. (p^) in deriving the identity (p7|). 

The Eq. (|25|) , though very simple, has the very attractive feature that it explicitly 
shows the contributions resulting from different types of dynamical singularities: the 
bound states, the resonances, the virtual states, and the left hand cut. The poles of 
the scattering amplitude contribute only to the /oca/ terms in the representation (p5|). 
The sum on the r.h.s. of Eq. ( pSf ) includes the poles on both the first and the second 
sheets. However, the poles on the sheet II are canceled by the corresponding zeros 
of the S-matrix on the sheet I, so that the representation (P^), which determines 
the scalar form factor on the sheet I, contains only the poles corresponding to the 
bound states. 
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(a) (b) 




Figure 3: The scalar form factor corresponding to the Jost function defined by 
Eq.(^): (a) a virtual state aX k = —iO.lfi, (b) a bound state at k = iO.lfi. The 
solid line is the complete result for |y47(ii^)|, the dashed lines show the result in 
the scattering length approximation. The particle mass is set to m = 0.5, so that 
E = k\ 



2.3 Examples 

To illustrate the above described formalism we consider a simple example of a non- 
relativistic scattering of a particle with mass m by the Yukawa potential V{r) = 
gexp {—fir)/r. The Jost function in the first order in the coupling constant has the 
form 

T{}z) = 1-Vln(l ) (28) 

tk /i 

For the sake of simplicity, we neglect the higher order terms because all analytical 
properties of the scattering matrix and the form factor are strictly preserved in any 
finite order of the Jost function expansion. [] The S-matrix corresponding to the 
Jost function (|28|) for g < has only one pole at /c = in: it is either a virtual state 
on the sheet II (/t < 0) or a bound state on the sheet I (/t > 0). The calculations 
using Eq. (|25|) are straightforward, and typical examples for the cases of virtual and 
bound states are shown in Fig. ^ 

One commonly used approximation is to neglect explicit effects of the left hand 
cut and consider only a few poles of the scattering amplitude close to the physical 



"'^The expansion of J-{k) in powers of the coupling constant is straightforward. 
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region. In the scattering length approximation the S-matrix has the form 



where a is the scattering length. In this case we have only the pole term in 
at 2mE = k"^ = a~^, and the subtraction term vanishes for A — oo (^(oo) = — 1). 
This gives the well known result: 

^'^^^^ = 1 (30) 
1 + ik[h)a 

A comparison of this result with the above considered example using a = is 
shown in Fig. ^. The range of validity of the scattering length approximation for 
a weakly bound state is limited by the energy range \E\ ~ 2m/a^. Note, that the 
approximation Eq. (^) does not have the correct asymptotic behaviour: A(E) 
1. 

Another useful example is the pole approximation for K matrix. In the single- 
pole approximation, the K matrix has the form 

K{E) = (31) 

and the corresponding S matrix is 

_ l+^kiE)KiE) 

- l^Jk(E)K{E) ■ ^ ^ 

Assuming Ei > and < 71 < 2^/Ei, one finds two resonant poles at E = Er and 
E = E; where Er = (-^72/2 + ^/K^^l[/4:Y . Using Eq. (|2|) we get 

{E,-E + zkj{E)^,) C{E-E,) 

2{E,-E) {E-E,){E~E;) ^ ^ 

where a is the residue of the form factor at E = E^ and C is a normalization 
constant. At 71 7^ 0, the form factor has no pole at the bare pole position E = Ei, 
therefore Ez = Ei, and the form factor is given by 

ME) = - E.:)(E - e;) '^'^^ 

where we set C = 2Ei. The r.h.s. of (^) has no poles on the first energy sheet 
because the nominator {Ei — E + iki{E)'yi) = {Ei — E — ikjj{E)'-^i) vanishes at 
E = Er,E*. It is easy to verify that the result (R3) is equivalent the standard one 



pole approximation |27 



ME) - (^._^!'.,,(^),.) (36) 

This result can be easily generalized to the case when K matrix has more than one 
pole. 
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3 The analytic form factor and the S matrix in 
the complex s plane 

In the previous section we have introduced the new representation for the form 
factor in the complex k-plane in the language of scattering theory. It is of course 
straightforward to re-express all the results in the above section in relativistic quan- 
tum field theory, in the complex s plane (here s denotes the center of mass energy 
squared). Hence the spectral representation of the form factor can be written down 
using the well known LSZ reduction formalism, 

ImA = ApT* , (37) 

where p denotes the kinematic phase-space factor. Taking vrvr scattering for example. 



p=^l-Aml/s. (38) 
In Eq. (|37|) T denotes the scattering T matrix which satisfies the optical theorem, 

ImT = TpT* , (39) 



for the physical value of s. Discussions parallel to that in section p.2| can be made 
and Eq. ( |25| ) can be recasted as. 



Aj{s')disc 



l+Sijs') 
Si{s') 



where the position of zeros of 5*/ on the complex s plane are denoted as Zi, and 
Sj are the position of bound state poles of Sj (and A). The P„(s) is a (n-l)th 
order real polynomial when A{s) obeys a n-th subtracted dispersion relation. The 
discontinuity of S on the l.h.c. L = (— oo, 0] manifests itself in the left-hand integral 
on the r.h.s. of the above equation, even though A itself does not contain left-hand 
singularities. Possible (n-th) subtractions on the integral in Eq. ( ^0]) are understood. 

The Eq. (^0]) offers a convenient and powerful expression for studying the analytic 
structure of the scattering 5* matrix and form factors. For example, it sets up a 
relation between the form factor in the time-like region and in the space-like region, 
provided that the S matrix is known. However, we will leave it aside (except for 
the analysis made under some simple approximations discussed in sec. |2.3D and 
go directly to discuss the analytic property of the scattering T matrix. In fact, 
it is easy to understand that the T matrix itself satisfies Eq. (^) with only minor 
modifications. That is, unlike the form factor, T also contains left-hand singularities. 
After some algebraic manipulation one can obtain (for the detail of the proof, see 
Appendix) , 

= t4^{E , , V + E ^ + Pniz) 

1 + S ^ 2p{Zi)S'{Zi){z - Zi) ^z-Sj 

d^sc frr^'^i±^ 

4 ' ^ 



2m 
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The relation Ti{zi) = i/{2p{zi)) is used in deriving Eq. (P|). If we define F = 
T(l + S')/5' so tlie above equation can be rewritten as, 



The Eq. ( ^4 ) sets up a dispersion relation for F which contains no right-hand cut. 
We can further express the S matrix in terms of F, 



Apparently F is real when z real and greater than 4m^. Actually ReT = -j in 
the physical region, hence F is the analytic continuation of 2ReT on the entire cut 
plane, and now it becomes obvious why in Eq. (|^) there is no right-hand cut. Also 
we have, 



Im^T = ^ , (44) 

for physical value of s. It is worth pointing out that Eq. (|^) relates the imaginary 
part of a given partial-wave amplitude on the unitarity cut to the imaginary part 
of the same partial- wave amplitude on the l.h.c, in a complicated manner. It 
is interesting to compare Eq. (p^ with the following relation obtained using the 
Froissart-Gribov representation for the partial wave projection. 



1 -|_ (^lY+L 

Im,T/(.) = ^-^-^EE(2^' + 1)^S? 



X 



s — Ami 



which p8[ relates the left-hand discontinuity to the right-hand discontinuity in a 
linear way, but all partial-waves are involved. To further analyze Eq. (P5|), we notice 
that \pF{s)\ < 1 in the physical region, since in the physical region 

pF{s) = sin(25) . (46) 

It is therefore easy to understand that at values of s when 6 = 7r/4 + n7r/2 we have, 

^(PF)=0, (47) 

and 

^ipF) < (pF = 1) ; ^(pF) > (pF = -1) . (48) 

These requirements lead to additional constraint on the parameters in the expression 
of the S matrix. Especially Eq. (|47| ) and Eq. (^Sf) indicate that the pole parameters 
and the discontinuity on the left are correlated to each other. For example, without 
further knowledge on the l.h.c. integrals, Eq. (|47|) (Eq. (|48|)) enables us to establish 
a sum rule (inequalities) that the left-hand integral has to obey. 
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The Eq. (^) is important in the sense that it relates the experimental observable 
(the partial-wave phase shift) to the pole contributions and the l.h.c. contributions 
in a simple and elegant manner. Unlike the conventional K matrix approach, the 
pole parameters are physical and the two contributions from poles and l.h.c. are 
additive. As we will see in the next section that we will be able to calculate the 
l.h.c. integral in Eq. in a reliable way, hence Eq. ( ^^D offers a convenient method 
to fit the pole positions. It is worth emphasizing that Eqs. (^21) and (|46|) are valid 
for any partial-wave amplitude hence they afford a unified approach for determining 
resonances in different channels. This is remarkable since the output from different 
channels can be compared to each other and this feature is very important to evaluate 
the quality of the method and approximations being used in the fit. 



4 Estimates on the left— hand cut effects of nn 
scatterings in chiral perturbation theory 

In above discussions we have re-expressed the form factor and the S matrix in 
a way that their dependence on isolated singularities and the l.h.c. are explicitly 
exhibited. Since the effects of the unitarity cut are dissolved, the impact of the left- 
hand singularities on the analytic structure of the amplitude is singled out. These 
formulas can be particularly useful when the fine structure of the l.h.c. becomes 
important, which has been ignored by most applications of the K-matrix approach. 
For example, when determining the pole position of the a resonance from tttt phase 
shift, since the a particle's mass is rather low, one has to carefully take into account 
the fine structure of the l.h.c, as being emphasized recently in Ref. |T^, Inspired 
by this we in the following carefully discuss the effects of the left-hand singularities. 
In order to achieve this and to reduce the model dependence as much as possible 
we need a method to estimate the l.h.c. of nir scattering amplitudes in a reliable 
way. The chiral perturbation theory (xPT) affords us such a method. The xPT 
treats the T matrix as a perturbative power expansion of the external momentum 
and works very well at low energies. Of course, a calculation in xPT to any finite 
order in the power expansion would not reconstruct the pole structure, therefore 
xPT encounter problems when there is a resonance near the physical region on the 
right-hand side of the complex plane. However we argue that its predictive power 
in the region near the l.h.c. would suffer much less from such a problem since the 
pole position is further away from the left-hand cut. Our strategy is clear from the 



above discussion: We extract from the 1-loop xPT results of the T matrix |^0[ the 
term relevant to the left-hand discontinuity to estimate the left-hand integral in 
Eq. (^). To be precise, the 1-loop xPT results for tttt scatterings are, 

T^=o = 3A{s,t,u) + A{t,u,s) + A{u,s,t) , 
T^=i = A{t,u,s) - A{u,t,s) , 

T^=2 = A{t,u,s) + A{u,t,s) (49) 
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where 

A{s, t, u) 
Bis, t, u) 



s — m 



^ + B{s, t, u) + C{s, t, u) + 0{E^ 



P 

J TT 

^{3(^2 - mt) J(s) + [tit -u)- 2mlt + Amlu - 2mi\ J{t) 
+ {t^u)} , 



(50) 



and the function C is a polynomial of s, t and u which is irrelevant here. The 
function J(s) is defined as, 



J{s) 



1 



167r2 



pin 



+ 2 



(51) 



The partial wave projection can be carried out: 

1 /-o 



T'j{s) 



327r(s - Ami) - 

/4m| 



/ 

Jin 



dtPj(l + 



2t 



Ami 



)T'(s,t,u) 



(52) 



from which the discontinuity of T on the left can be obtained under certain condi- 
tions [EHl, 



ImiTi(.) 



1 + (-1)^+*^ f*"'' 
327r(s-4ml) J Ami 



dtPj(l + 



2t 



s — Ami 



)ImT/(s,t) ;s < . (53) 



However, for scattering amplitudes in chiral perturbation theory, the left-hand cut 
can be directly extracted from the analytic expressions of partial wave amplitudes. 
In the I = J = channel the result is, 



im^r°(.) 



1 



15367r2/4(s - Ami 



-{2 -In 



Ami 



Ami — s + 



(25m^ - 6m^s) 



+ ^-s{Aml-s)(^-s' - ^mls + 25mt)} . 



(54) 



However, this result is still not directly applicable here. Remember that the function 
F is an analytic continuation of (2 times) ReT, which is the real part of T defined 
in the physical region. Since Im^^T also develops a discontinuity on the left this part 
must be subtracted from Im^T. Therefore the correct expression for Im^F is. 



ImLF = 2(ImiT - Rcilm^T) 



(55) 



The 1-loop result of Im/jTQ is simply obtainable using the Born term amplitude and 
the optical theorem. 



Im^To = p 



'2s — ml 



(56) 



for which we have Re^Im/jT = Im^T since p is real on the left. It is worth pointing 
out that in Eq. (|55|) the main contribution to Im^^F comes from the second term on 
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the r.h.s., f\ i.e., the term Im^T is numerically rather small comparing with Im^^T 
evaluated on the l.h.c. (For previous discussions on the smallness of Im^T, see for 
example Ref. fl^.). 

We also list the results in J = J = 1 and / = 2, J = channel from chiral 
perturbation theory: 

-s(4m^ ~ -5) + 12m^(6m^ + 13m^s — 3s^) 



xln ^, : , (57) 

'Ami — s + V— s 



967r/2 



+2{ml + 3mts) x In ^ " } , (59) 

/ 4j7^2 _ s -|- 




ImnTi = p\——f \ , (60) 



where /^r = 93.3MeV and the expressions listed above are valid up to 0{p^) term 
in the chiral expansion. It is worth pointing out that at s = 4m^ and 2m^ 

the lowest order partial wave amplitudes Tq , Tl and Tq vanish respectively, as a 
consequence of the Adler zero condition. 

5 The determination of the a resonance in tttt 
scatterings 

5.1 A description to the fit procedure 

Having obtained the analytic expressions of Im^F for the partial wave amplitudes 
it becomes possible to use the experimental data on tttt scattering phase shifts to 
determine the pole positions in various channels. We fit both IJ=00, 11 and 20 
channels. In the following we briefly describe the procedure and the method for the 
fit: 

1. We assume the tttt scattering amplitude T satisfies a once subtracted disper- 
sion relation from physical considerations. So our left-hand integral is once 



^The necessity to have ImflT included in Eq. ( p5[ ) can be evidenced by Fig. |^: without the 
contribution from Ini^T the slope of line A or B from xPT would be significantly smaller in 
magnitude. See later text for more discussions on Fig. ^ 
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subtracted at s = s = 4m^ and s = 2m^ for IJ=00, 11 and 20, re- 

spectively. The positions of the subtraction points are chosen only for the 
convenience of discussions. 

2. A dispersion integral in xPT needs many subtractions and this is artificial 
because of the bad high energy behavior of xPT amplitudes. So we always 
truncate the once subtracted integral at — A^p^ where A^py ranges from, for 
example, 700MeV to IGeV. We argue that the influence from the region \ y/s\ > 
A-^PT is negligible to the physics we are concerning. We realize that the value 
|s| ~ IGeV is too large for chiral perturbation theory to be valid. We however 
take this value just for the purpose of testing to what extent our fit results 
depending on the behavior of ImiF. We have to make sure that the mass and 
the width of the resonance obtained under such an approximation have to be 
insensitive to the choice of the cutoff. To compare with perturbation results 
we will also extract Im^F from the unitarized T matrix, i.e., the [1,1] Pade 
approximant, in the I=J=0 channel. 

3. In the three channels IJ=00, 11 and 20, we assume one pole in the first two 
cases. In the IJ=11 channel the threshold behavior (T ~ near threshold) 
is considered to reduce one parameter (the subtraction constant). In our fit 
each pole contains 4 parameters. Two of them are related to the position of 
the pole. The other two are related to the couplings of the pole to F. That is 
the coefficient aj = i/ {2p{zi)S'{zi)) in Eq. (|42|) . One may relate to Zi in the 
narrow resonance approximation or in models. But for the wide resonance a 
there is no simple relation between the two, so we take aj as completely free 
as an approximation. However, we also fit in the I=J=1 channel by treating 
ReOj and/or Imoj as the one obtained from the narrow width approximation. 



5.2 Numerical results and discussions 

The numerical results of the fit in the IJ=11 and 20 channels are presented below: 



1. The I=J=1 channel: The data are taken from table VI of Ref. [^|. We 



find that the l.h.c. effect is very small in this channel. When taking ap for 
free (so it is a 4 parameter fit) we obtain Mp = 753MeV, Tp = 142MeV 
for A^PT = l.OGeV. Without l.h.c. while taking ap for free, we have instead 
Mp = 756MeV, Tp = 144MeV (See Fig. | for the fit). We also tried to fix Rea^ 
and/or Ima^ using the narrow resonance approximation. The pole position 
obtained using different approaches agree with each other within a few MeV. 
The is sensitive to different approaches, but the pole position is not. 

2. The 1=2, J=0 channel: 

Since there is no resonance pole in this channel it is a one parameter (the 
subtraction constant) fit here when A-^pt is held fixed. We notice that the 1=2, 
J=0 channel favors a large value of A^pt- As we see from the line B of Fig. ^ 
the phase shift is reproduced ideally by the l.h.c. integral plus one subtraction 
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constant. The scattering length is estimated to be Og = —0.052 which is about 
2(T away from the experimental value, al{exp) = —0.028 ± 0.012 f\ and the 
most recent value from the E865 Collaboration |35|: a^lexp) = —0.036 ±0.009 
The case without the subtraction constant is also depicted in Fig. ^ as line 
A. The latter is equivalent to imposing the Adler zero condition for the 1=2, 
J=0 partial-wave amplitude, since in here the dispersion integral is subtracted 
at s = 2m^. The scattering length Oq = —0.032 for the latter case, which is 
within la error bar comparing with the experimental value. See also footnote |^ 
for further information. 

In the I=J=1 channel the p resonance almost saturates the experimental phase 
shift and hence it demonstrates that the contribution from the left-hand integral 
must be very small, and it is correctly predicted by chiral perturbation theory. In 
contrast, no resonance exists in the 1=2, J=0 channel. It is of course very satisfiable 
to see that the contribution from the left-hand integral to 6q is in the right direction, 
and combining the contribution from the subtraction constant it can saturate the 
phase shift. 

From the above discussions we see that gives satisfactory results in both 

the IJ=11 and 20 channel. We may obtain an impression from the above discussion 
that the chiral prediction works well at least in qualitative sense, i.e., the order of 
magnitude and the sign of the left hand integrals. In the following we turn to discuss 
the more interesting case of the I=J=0 channel. 

In an earlier version of the present paper, we take the data the same as that 



in Ref. [0], that is the CERN-Munich data [32, 33] combining with the data from 



i^e4 decay near threshold 153] obtained more than 20 years ago. The most recent 



experiment on Ke4: decay performed by the E865 Collaboration ]]35[ is remarkable. 
It affords high statistic data on the vrvr scattering phase shifts near threshold from 
which the scattering length parameter Oq can be determined with much improved 
accuracy: Oq = 0.228±0.012±0.003. We therefore incorporate the new experimental 
results in our fit below. Here we truncate the CERN-Munich data at 900 MeV in 
an attempt to reduce the pollution from the /o(980) resonance. Since the latter 
(or more cautiously, its second sheet pole) is only a narrow resonance, its influence 
should not be very important to the determination of the a resonance. As we see 
from table 1 the pole position is not very sensitive to A^py. Including the l.h.c. 
contribution does reduce the mass and width of the a particle, but the effect is not 
dramatic. The ^ire good for all cases, see Fig. |^ for the fit. It is very impressive 
to notice that in the I=J=0 channel the contribution from the left-hand integral 
to sin(25Q) has the wrong sign comparing with the experimental value. To clearly 
show this we draw in Fig. ^ the contribution from the left-hand integral alone to 
sin(2(5Q) in several cases, the curves are always negative and concave regardless of 
the different choice of cutoff parameters. The procedure presented above clearly 

■^AU the experimental value and the results from chiral perturbation theory are taken from 
Ref. unless quoted explicitly. 

^This is in qualitative agreement with the result of Ishida et. al. | |3^ from a linear a model 
calculation. 
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J- a 


A^PT = 1000 


477 


564 


A^PT = 850 


503 


587 


A^PT = 700 


519 


605 


Fade 


519 


579 


no l.h.c. 


544 


607 



Table 1: The pole position of the a resonance obtained by a fit of the data from 
], ^ and . The scattering length Oq obtained in the global fit ranges 
0.26. The data in the table are in units of MeV. 



Refs. P 



from 0.25 



demonstrates the existence of the a meson, in a model independent way Q since 
within our scheme the 'background phase' or the cut becomes calculable, from the 
principle of maximal analyticity no other contribution except the a pole can be 
introduced to fit the experimental data. The subtraction constant, which represents 
the contribution from high energies, can not generate the convex curvature that the 
experimental value of sin(2(57r) exhibits. 

A cautious reader may worry about the convergence problem of the chiral ex- 
pansion in the large A^p^ region, even though the qualitative behavior of the left 
hand integrals obtained from perturbation theory can be examined by varying the 
cutoff parameter. To overcome the difficulty we also estimated IovlF in a 'uni- 
tarized approach', though crossing symmetry is no longer maintained in such an 
approximation. That is we use the [1,1] Fade approximation to obtain the unita- 
rized scattering amplitude (in which the standard values of the Lj parameters of the 
Gasser-Leutwyler Lagrangian are used) and extract the quantity Im^^F on the left, 
and insert the obtained Im^F into Eq. (^2]). The left hand integral evaluated from 
the unitarized amplitude is also plotted in Fig. ^ and we find that no major conclu- 
sion is changed. It is worth pointing out that such a use of the Fade approximation 
is different from the standard usage of unitarization (see for example Ref. |p0|-|2^). 
Essentially we only need the information of the unitarized amplitude on the cuts 
evaluated on the left, and the two approaches are not technically equivalent. 

However, the scattering length parameter Og obtained in the global fit given 
in table 1 ranges from 0.25 ~ 0.26, which is about (or more than) 2a away from 
the newly obtained experimental result (It may be worth noticing that without the 
new data the Og parameter obtained in our procedure would be much larger, 
ag = 0.35 ~ 0.37.). To remedy this we also include in our fit program by brute force 
the experimental constraint Og = 0.228 ± 0.012 ± 0.003. The results are given in 
table 2. The fit result for Og, ranging from 0.230 ~ 0.231, is stable against different 
treatment of the left-hand integrals. 

By comparing table 1 with table 2 we find that the pole position of the a res- 
onance is sensitive to the scattering length parameter, but the infiuence is not as 
dramatic as in the case without the new KeA data from the E865 Collaboration: the 



^ To be more careful, one may add: provided that a chiral perturbation theory estimate on the 
l.h.c. is qualitatively correct. 
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A/f 


J- cr 


A^PT = 1000 


454 


648 


A^PT = 850 


479 


658 


A^PT = 700 


500 


662 


Fade 


498 


642 


no l.h.c. 


525 


654 



Table 2: The pole position of the a resonance with the new experimental constraint 
on Oq. The value of Oq in the global fit ranges from 0.230 ~ 0.231. 



pole position would be much more sensitive to the scattering length parameter in 
the latter case. 

In discussions above we limit ourselves to a single channel analysis, therefore 
the possible influence from higher resonances, especially the /o(980), is omitted. 
We point out here that within the single channel approximation the influence of 
the /o(980) resonance to the fit can also be estimated (to be exact, in a model 
dependent way) using the method of Ref. [^], by subtracting the phase contributed 
from /o(980) alone. The mass and the width of the a particle change only slightly 
when /o(980) is taken into account, but the qualitative role of the l.h.c. is still 
unchanged, that is it has only a mild influence. 

As we can see from tables 1 and 2 that the pole position of the a resonance on 
the complex s plane, s^, moves towards left once the value of aj] is decreased ^ The 
phenomenon that the real part of the pole position of the a resonance is small has 



been noticed and investigated by Anisovich and Nikonov ||T8[. An interesting mech- 
anism was proposed to explain this phenomenon by suggesting a strong singularity 
associated with the l.h.c, which is simulated by a series of poles on the negative 
real axis fixed by the N/D equation and the experimental data. Our calculation 
confirms that the inclusion of the l.h.c. effect pushes the a pole towards left but the 
effect is mild, and Re[so-] is also found to be sensitive to the value of the scattering 
length parameter. The difference between our treatment on the l.h.c. and that of 
Anisovich and Nikonov is that the left-hand singularity of the T matrix in our case 
comes solely from the absorptive singularity, i.e., the 27r cuts in crossed channels. 
Due to relativistic kinematics, the quantity F contains an additional contribution 
from the s channel absorptive singularity, which is also from the 2tt cut. [] It may 



be illuminating to quote from Ref. [^|, "There are no such singularities which are 



associated with the quark-gluon structure of hadrons, since there are no absorptive 
thresholds related to this structure" . 



^In the previous version without the newest K^.^ data, the fit value of Re[S(j] is much more 
sensitive to the value of ag. Even a negative Re[so-] can be obtained when ao(x-Pr) is used in the 
fit. 

^ Related discussions on the left-hand cuts of tttt scattering amplitude in the recent literature 
may be found in Refs. |^ and 
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6 Conclusion 



The method proposed in this paper to discuss vrvr interactions is based on a dispersion 
relation set up for the analytic continuation of the real part of the scattering T 
matrix. The main formula Eq. (|42|), though very simple, is shown to be very useful 
in clarifying different contributions from poles or the left-hand cut to the scattering 
phase shift. The procedure in deriving Eq. (^) fully respects analyticity, unitarity 
and crossing symmetry. 

Our estimate on the left-hand cut of the scattering T matrix obeys the standard 
rule of S matrix theory, that is the left-hand singularity comes solely from physical 
absorptive singularities in the crossed channel. In here they are the 2n cuts from t 
and u channels. For the quantity F, it also contains the part originated from the s 
channel absorptive singularity, as an effect of the relativistic kinematics. 

We have carefully examined the reliability in using chiral perturbation theory to 
estimate the left-hand cut effects in various channels. Chiral perturbation theory 
encounters the problem of a bad high energy behavior, especially in the I=J=0 
channel. This drawback was partially corrected by using a 'unitarized approach' to 
the perturbation series. It is remarkable to notice that the l.h.c. contribution to 
sin(2(5Q) is negative and concave which clearly demonstrates the existence of the a 
resonance, according to the principle of maximal analyticity. 

We have shown that the left-hand cut effects are mild in determining the pole 
position of the a resonance, and the effect of the scattering length in the determina- 
tion of the a pole position is also clarified. The estimated central value of the mass 
and the width of the a meson are found to be 498 MeV and 642 MeV, respectively, 
if we are allowed to quote from the Fade solution with the experimental constraint 
on Qq. The new experimental data from the E865 Collaboration on KeA decay is 
found to be crucial in minimizing the uncertainty caused by aj] in the determination 
of the a pole position within our scheme. 

The current procedure can be extended to discuss the more general case of 
coupled-channel system where the /o(980) resonance has to be taken into account |^ 



No major conclusion on the pole position of the a resonance and the role of the left- 
hand cut and the scattering length parameter is changed. 
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A Appendix 

By definition 

^ = 1 + 2ipT (61) 
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and from unitarity relation 

ImT = TpT* , (62) 

We can get the analytic continuation of T on the second sheet by using the reflection 
property, 

T^{s + ie) - T^\s + ie) = 2ip{s + ie)T^(s + ie)T^^(s + ie) , (63) 
or more concisely, 

rpl rpl 

T'' = - = — (64) 

Suppose the integration over the infinite circle is zero and using Cauchy's theorem 
we get 

2m Jcr s' — s 2m Jcl s' — s j s — Sj 



TT J4m| S' — S 2m J-oo s' ~ S ^ S — Sj 

TT J 4ml S' — S 2m J-oo s' — S j S — Sj 



On the r.h.s. we have 



p{s)T\s)T'*{s) = p{s + ie)T'{s + i€)T'\s + te), 

= -p{s - ie)T^\s - ie)T\s - ie) , (66) 



therefore we have 



^ J_r p{s')T{s')T^\s') ^ 1 |0 disdm^^,^^ 
277 Jcr s' — s 2m J-oo s' — s j ^ ~ 

1 f p{s')T{s')T''is') , , If p{s')T{s')T''{s') ^ , 



3 '-'j 



r p^s')ns'Jl-^s')^,_^l 
2tt Jc s' — s 27r Jcl s' — s 



-ds' 



J 



where 



2m J-oo s' — s ^ s — Sj 

2m Jc s' — s 2m J-oo s' — s j s — sj 

t + ^r fec(T(.')^)_^^,^^ (67) 

27ri J-oo 5^ — 5 ~^ S — Sj 

3 



-i Jc 



2m Jc s' — s 



2S^ 2S^{s,)'{s, -s) \ 2{s - Sj) ■ ^^^^ 
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In the derivation of Eq. (|B7D, by combining the second and the third term on the 
right hand side of the second equation of (^TD and using 

we get the third equation of (^). And by substituting $ into ( 167D we obtain 

S 



T s) 



S 



Y 2p{s,)S{s,)'{s -s{) y s- s, 2m I - ^' - ^ 



■oo S — S 



(70) 
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Figure 5: Contributions to the scattering phase from the left-hand integral (once 
subtraction at s = 2m^ is used) in the / = 2, J = channel. Line A corresponds to 
setting Tq (s = 2m^) = and be parameter free, and line B corresponds to a free 
subtraction constant fixed by minimizing . The data are from Ref. ||32|| . 
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Figure 6: A typical fit of 5 parameters (4 resonance + 1 subtraction constant) in tlie 
I=J=0 cliannel, witli A^pT=0.7GeV and the experimental constraint on Oq. Different 
estimates on the left-hand integral are also plotted: The dotted line corresponds to 
^xPT — O-QGeV, the dashed line corresponds to A^pt = 0.7GeV, the dot-dashed 
line corresponds to A-^pt — 0.85GeV^ and the solid line corresponds to the Fade 
solution. 
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